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1. Introduction 

1.1 Chiral Perturbation Theory on the Lattice 

One of the most important uses of Chiral Perturbation Theory (#PT) is the extrapolation of 
results from Lattice QCD simulations. Finite volume effects, discretization errors and quark mass 
extrapolations can all be addressed using #PT. In this context, #PT can be regulated with either a 
continuum formalism or a lattice formalism Though the lattice approach typically leads to 
more complicated algebra, it has the advantage of being directly amenable to numerical methods. 
A recent study of volume dependences in 1-loop lattice #PT demonstrates how simple numerical 
methods are able to replace extensive algebra Besides its relevance to QCD simulations, lattice 
^PT is vital for nuclear simulations with chiral effective theories 

In lattice regularized #PT , the meson fields live on the lattice sites, and gauge fields as usual 
on the links. Derivatives in the Lagrangian are replaced by appropriately covariantised finite dif- 
ferences. Lattice Chiral Perturbation Theory has been used to compute baryon magnetic moments 
[Q] and to simulate neutron matter [Q]. 

Since Feynman rules become much more complicated on the lattice due to the appearance of 
trigonometric functions and the loss of Lorentz symmetry, deriving them manually becomes an 
extremely tedious task, making an automatic procedure desirable. 

1.2 Automated Generation of Feynman Rules 

Procedures to automatically derive Feynman rules for lattice gauge theories from their lattice 
actions have been known since [§]. A new, efficient and adaptable method was presented in [0]. 
The great flexibility of this latter method makes it possible to adapt it to the case of #PT on a lattice 
without much difficulty by adding support for separate left- and right-vector fields and extending 
the SU(A0 algebra to include quark mass matrices. The structure of a meson «-point vertex with 
mesons of flavours a, incoming with momenta pt is 



where each term can be described by an entity (j,X, T,x,y, {vj}), with F the usual #PT parameter 
(pion decay constant in chiral limit), % the chiral order and / the amplitude of the term, T a SU(A0 
trace, and x, y and Vj lattice sites. Entities can be multiplied with the multiplication rule 



following from the above representation of the vertex. 

Fields can be expressed in terms of basic entities, and arithmetic operations on fields can be 
turned into operations on the entities in their expansions. Once the field algebra has been imple- 
mented, the perturbative expansion of the action expressed in terms of the fields into entities will 
occur automatically. 

The expansion algorithm is wrapped into a PYTHON program called CHIRPY, which performs 
the expansion of an arbitrary lattice #PT action given as input and produces the Feynman rules 




(f,X,T,x,y,{vj})*(f',X',T',x',y , ,{ V / f }) = 
(ff',X + X',T* r,x,y + y'-x', { Vj } U {v' y +y-x'}) 
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in a human-readable LTpjX, compileable Fortran or special machine-readable format (Fig. The 
relevant flavour, Lorentz and Taylor |§] algebra is provided in a companion Fortran 95 library called 
CHIRPER. Other output formats (e.g. for analytical evaluation by a computer algebra program) 
could be added easily as required. 

1.3 Automated Generation of Feynman Diagrams 

There exist a number of well-known packages for the generation of continuum Feynman dia- 



grams, such as FeynArts Q9p or Qgraf plOp. None of these, however, easily support the vertices and 
counterterms of arbitrary order that appear in lattice #PT. 

For this reason, a straightforward implementation of Wick's theorem is implemented as a 
Python program called Madeline, which creates the Feynman diagrams contributing to a spec- 
ified «-point function at a specified chiral order. The diagrams are output both as LATpX/FeynMF 
figures and as Fortran functions using the CHIRPY/CHIRPER routines (Fig. ||). 

Using the integration routines provided by the CHIRPER library, these functions can be in- 
tegrated numerically for an almost automated evaluation of the corresponding lattice Feynman 
diagrams (where internal propagators carry external on-shell momenta, some user intervention will 
be necessary to assure the choice of the correct integration contour in the complex plane). 



2. Usage Overview 



# of legs and loops 



F90 Feynman rules 



CHIRPER 
library 



.F90 diagram code 



LaTeX/FeynMF 



ChiPT application 



Figure 1: Schematic of usage of the tools 

Fig. U shows a flow schematic of our tools: The user specifies the action as an input to the 
CHIRPY main script, which outputs the Feynman rules in the format chosen by the user. 

To compute a particular quantity, the user specifies the number of legs and the loop order as 
input to the MADELINE script, which outputs representations of the relevant Feynman diagrams in 
both Fortran code and as graphics (where the LTpX/FeynMF output may have to be edited by hand 
to achieve an esthetically pleasing result). 
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, „ TTTT ,™, ,17 j t „„ a o rr, nnnr MODULE XPYJ/ERTEXJ14 

1 CHIRPY run on Wed Jun 22 16:48:57 2005 IMPLICIT NONE 

«*(*«) ,-(*«) vertex: CONTAINS 

FUNCTION VERTEX_M4 (k , f , chi.Order) 



iL ira (*« + 4 2 >+*<'>-t< 2 >) 



USE chp_generic 



2 v ' USE chp_momenta 

+L 1 4cos(2t< 1 >)} USE chp.taylors 

H (*«) (SO) *- (*<«>) vertex: IMPLICIT NONE 

TYPE (mom) , INTENTCIN) :: kC4) 
£ + + INTEGER, INTENTCIN) :: f(4), chi.order 

TYPE(taylor) : : vertex_m4 
INTEGER :: flavkey, i,niu,nu,rho,sigma 



1 cos (kf) - 1 cos fe>) f lavkey = 

2 / ,,\ 2 / D0 1=1,4 

jL4iJ cos (2t<»>J- -£44003(2*1?) J flavkey = flavkey + f Ci)*4**C4-i) 



ENDDO 

rau=0 ; nu=l ; rho=2 ; sigma=3 



sfc 1 ' + *< 2 > + kf) SELECT CASECchi.order) 

CASEC2) 



SELECT CASECflavkey) 



--L.M? cos (V 2 >-*< 3 >1 CASEC90) 

2 ' v"» " ) vertex_m4 = k ! pi+CkCD) pi+(k(2)) pi-(k(3)) pi-(k(4)) 

--L,iia*s(-k<»+kf) 16.D0/3.D0+ 2. DO/3. DO * mass_pip**2 k 

l r , / , m M \ , , /„,,„ „ M \ - 4. DO/3. DO * TCOSC (/ k(l) /),C/mu/)) k 

-- £ . 1 4cos(- t « +t P)) +£4lira (2*« + 2 1 f)) . 4 M/3 M . TC0S( u k(1) /)](/nu/ „ & 

+L4i;cos(2*< 2 )+2tj?>) - 4. DO/3. DO * TCOSC (/ k(2) /),(/mu/)) k 

WoJ / 2t (i, + t e) + ,.< 3 A " 4.D0/3.D0 * TCOSC (/ k(2) /),(/nu/)) k 

-L4i;cos(*m + tji 2 > + 2*( 3 >) - 4. DO/3. DO * TCOSC C/ kC3) /),C/nu/)) & 

+ te + 2U£\ cai (>f» + m - 4 - D °/3- D * TCOSC C/ kCl) /),C/rho/)) & 

+ l + 4 *' ^» + ' ) - 4.D0/3.D0 * TCOSC C/ kC2) /),(/*./)) Jt 

+ (2 + 2L4iJ)cos(4 2 > + 4 3 )) - 4. DO/3. DO * TCOSC C/ kC3) /),C/rho/)) 



■ 4. DO/3. DO * TCOSC C/ kCl) /), C/sigma/)) Jt 

■ 4. DO/3. DO * TCOSC C/ kC2) /), C/sigma/)) & 
-ii 1 cos(*< 1 ) + t< 2 >-t< 1 >+t< 2 >') - 4. DO/3. DO * TCOSC C/ kC3) /), C/sigma/)) & 

1 ; : + 2 * TCOSC U kCl) + kC3) /),C/mu/)) k 

cos (i< 2 > + i< 3 > + *») - *WJ + 2 t TCOSC C/ kCl) + kC3) /),C/nu/)) It 

+ 2 * TCOSC C/ kC2) + kC3) /),C/mu/)) k 

Figure 2: (left) I^TgX output generated by CHIRPY in human-readable mode, (right) Fortran output gener- 
ated by CHIRPY in compileable mode. Shown are only the first pages of the multi-page outputs. 



The generated code interfaces with the CHIRPER library and the Feynman rules generated 
by CHIRPY in either compileable or machine-readable format. With just a small amount of user- 
written code to bind these components together, a fully functional application of lattice #PT can be 
produced. Easy-to-use interfaces to initialisation and integration routines are provided by CHIRPER, 
making the user-written code rather straightforward. 

The structure of the interface between CHIRPY and MADELINE allows to reuse the same 
Feynman rules for a different set of diagrams and vice versa. 



3. An Application to Finite- Volume Effects 

As an example of an application of this framework, we have repeated the one-loop computation 
of the volume dependence of the pion mass from [|J]. The CHIRPER software is so structured that 
it is easy to evaluate the Feynman diagrams using a variety of summation/integration methods, 



including VEGAS [11]. This allows to directly extract the difference between the finite- volume 
and infinite-volume values of the diagrams (instead of using a large, but finite volume to provide a 
value for the infinite-volume limit, as done in ||3|]). 

The plot of Fig. |] shows an example of the results which agree well with the earlier evaluation 
[Q] that was based on a hand-coded implementation of manually derived Feynman rules. 

The great reduction of the workload in terms of human time by the automated method makes 
a two-loop calculation appear viable, and a first direct calculation of the two-loop volume depen- 
dence of the pion mass is currently under way. 
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FUNCTION DIAGRAM_1 C k, p, f ) 
USE chp_flavours 
USE chp_momenta 
USE chp_taylors 
USE xpy_vertex_m2 
USE xpy_vertex_m4 
IMPLICIT NONE 
TYPE (mom) :: k(2) 
TYPE (mom) :: p(2) 
INTEGER :: f(2) 
TYPE(taylor) :: DIAGRAM. 1 
INTEGER :: hl,h2,h3 
DIAGRAM_1 = 0.D0 
DO hl=l,3 

h2 = anti(hl) 

DO h3=l,3 

DIAGRAM. 1 = DIAGRAM_1 + & 
1. DO/4. DO * & 

PROPAGATaRJCCltCl), hi) * & 
PR0PAGATQRJ12CCp(l)+p(2)-k(l)), h2) * & 
PR0PAGATaR_M2(k(2), h3) * & 

VERTEX_M4( (/ p(l) ,p(2) ,-k(l) ,-(p(l)+p(2)-k(l)) /) , Sc 
(/ f(l),f(2),anti(hl),anti(h2) /) , 2 ) * & 

VERTEX_M4( (/ k(l) , (p(l)+p(2)-k(l)) ,-k(2) ,k(2) /) , & 
(/ hl.l^antiaiS),!^ /) , 2 ) 

ENDDO ! h3 
ENDDO ! hi 
END FUNCTION DIAGRAM_1 



feCl) 

p(2) 





Figure 3: (left) Fortran code generated by Madeline, (right) FeynMF diagrams generated by Madeline. 



0.1 



\ 0.01 



Figure 4: Volume dependence of the pion mass at M n — 140 MeV 



4, Conclusions 

Chiral Perturbation Theory can be formulated on the lattice, and efficient methods from per- 
turbative Lattice QCD can be translated to the #PT case with only moderate effort. A software 
implementation of an efficient algorithm for generating perturbative expansions exists. The Feyn- 
man rules and diagrams can be generated in different formats suitable for use by either humans or 
computers. An almost automated computation of lattice Feynman diagrams in lattice regularised 
#PT is possible with these tools. 

The automated procedure has successfully been used to repeat finite-volume calculations pre- 
viously done manually. New finite-volume calculations at the two-loop level using the automated 
procedure are currently under way. 
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